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A  one-dimensional  model  of  activated  sludge  secondary  clarifiers  with  a  dispersion 

term  dependent  on  concentration  and  feed  velocity  was  developed.  The  model  provides 

predictions  of  effluent  and  underflow  suspended  solids  concentrations  and  sludge  blanket 

height.  Data  collected  from  a  full-scale  clarifier  at  the  Kanapaha  Water  Reclamation 

Facility  in  Gainesville,  FL,  were  used  to  evaluate  the  model.  Better  matches  to  observed 

concentration  profiles  were  achieved  with  the  current  model  than  with  a  gravity-flux- 

constraining  model.  In  addition,  the  model,  when  calibrated  using  concentration  profile 

data  from  experiments  in  which  the  sludge  blanket  reached  steady  levels,  successfully 

predicted  the  outcomes  of  the  five  experiments  during  the  test  period  which  exhibited 

continuously  rising  blankets.  These  failures  to  reach  steady  blanket  levels  were  not 

predicted  by  either  limiting  total  solids  flux  theory  or  the  gravity-flux-constraining  model. 


IX 


Since  concentration  profile  data  are  not  readily  available  under  normal  plant 
operating  conditions,  the  ability  of  the  model  to  be  calibrated  using  sludge  blanket  height 
data  was  investigated.  The  clarifier  model  was  coupled  with  an  algorithm  for  predicting 
sludge  blanket  height.  The  model  was  successfully  calibrated  using  blanket  heights, 
effluent  suspended  solids  and  return  activated  sludge  concentrations,  and  the  measured 
Vesilind  settling  equation  parameters  V0  and  b.  Model  validity  was  confirmed  by 
comparing  predictions  of  the  calibrated  model  against  separate  sets  of  data.  Out  of  forty 
clarifier  loading  tests  for  which  the  system  could  conclusively  be  determined  as  overloaded 
(a  clarifier  failure)  or  operating  acceptably  (a  successful  state),  the  blanket  height- 
calibrated  model  correctly  predicted  the  outcome  of  thirty-seven.  Two  experimental 
successes  were  predicted  as  clarifier  failures,  whereas  one  experimental  failure  was 
predicted  as  a  success.  In  contrast,  a  gravity-flux-constraining  model  and  limiting  total 
solids  flux  theory  predicted  success  for  all  cases,  and  thus  failed  to  predict  the  eight  cases 
of  clarifier  overloading. 


CHAPTER  1 
INTRODUCTION 


The  process  of  separating  solids  from  wastewater  effluent  in  the  secondary  clarifier  is 
critical  to  the  optimal  operation  of  activated  sludge  systems.  The  secondary  clarifier 
performs  two  functions  in  this  capacity.  It  thickens  the  sludge  to  a  high  concentration  for 
recycle  back  to  the  bioreactors,  and  it  clarifies  the  wastewater  effluent  reducing  effluent 
suspended  solids  and  effluent  biochemical  oxygen  demand  (BOD)  due  to  effluent  solids. 
Because  clarification  is  an  integral  part  of  the  activated  sludge  system,  it  is  useful  to  have  a 
reliable  clarifier  model  that  can  be  incorporated  with  existing  biological  process  models  for 
application  in  design,  optimization,  and  control  of  activated  sludge  systems. 

The  objective  of  the  current  research  was  to  develop  and  validate  a  clarifier  model 
capable  of  predicting  sludge  blanket  heights  and  effluent  and  underflow  suspended  solids 
concentrations  that  could  be  easily  combined  with  an  activated  sludge  biological  process 
model.  Full-scale  clarifier  loading  tests  were  conducted  at  the  Kanapaha  Water 
Reclamation  Facility  (KWRF)  in  Gainesville,  FL,  to  collect  data  for  model  development, 
calibration  and  validation.  The  settling  characteristics  of  the  KWRF  sludge  were 
determined  by  conducting  batch  sludge  settling  tests  in  parallel  with  the  loading  tests.  The 
experiments  were  conducted  during  three  intervals  over  a  six  month  period. 


The  developed  model  employs  a  dispersive  flux  term  in  which  the  dispersion 
coefficient  is  a  function  of  both  solids  concentration  and  influent  velocity.  The 
development  of  this  dispersion  function  and  calibration  and  testing  of  the  model  are  the 
topics  of  Chapter  2.  Initial  model  development  was  based  on  the  model  of  Hamilton  et  al. 
(1992)  which  has  a  dispersive  flux  term  with  a  constant  dispersion  coefficient  and  the 
model  of  Takacs  et  al.  (1991)  which  applies  a  constraint  on  the  gravity  flux.  A  simple 
implementation  of  the  model  with  a  dispersion  coefficient  dependent  only  upon 
concentration  was  calibrated  separately  to  three  sets  of  steady-state  concentration  profile 
data  reported  in  the  literature  (Pflanz  1969),  and  the  resulting  model  fits  were  compared 
to  those  of  the  Takacs  et  al.  (1991)  model.  The  model  was  also  modified  to  reflect 
structural  characteristics  of  the  KWRF  clarifiers  and  calibrated  to  nine  sets  of  steady- state 
concentration  profile  data  from  KWRF  clarifier  loading  tests  conducted  during  one  of  the 
experimental  periods.  In  analysis  of  the  model  fits,  it  was  found  that  model  performance 
could  be  substantially  improved  by  incorporating  a  dependence  on  influent  velocity  in  the 
dispersion  coefficient  function.  The  calibrated  model  with  the  concentration  and  feed 
velocity-dependent  dispersion  term  was  used  to  simulate  the  six  remaining  KWRF  clarifier 
loading  tests  from  the  same  experimental  period  (one  which  reached  steady  state  and  five 
which  did  not  reach  steady  state  due  to  clarifier  overloading).  The  model's  ability  to 
predict  clarifier  failure  due  to  solids  overloading  was  compared  to  that  of  both  the  gravity- 
flux-constraining  model  and  limiting  total  solids  flux  theory. 

The  requirement  for  concentration  profile  data  to  calibrate  a  clarifier  model  is 
problematic  since  these  data  are  not  generally  available  under  normal  plant  operating 


conditions.  Therefore  the  use  of  sludge  blanket  heights  and  effluent  and  underflow 
suspended  solids  concentrations  to  calibrate  the  model  was  investigated.  This  is  the  topic 
of  Chapter  3.  The  algorithm  for  calculation  of  sludge  blanket  height  from  Hamilton  et  al. 
(1992)  was  modified  to  yield  better  agreement  with  blanket  heights  observed  in  the  field, 
and  the  model  incorporating  the  modified  algorithm  was  calibrated  using  sludge  blanket 
height  data  from  nine  KWRF  steady-state  loading  tests  conducted  during  one  of  the 
experimental  periods  and  the  measured  sludge  settling  equation  parameters.  The  resulting 
model  parameters  and  the  experimentally  measured  sludge  settling  equation  parameters 
were  applied  in  simulations  of  steady- state  loading  tests  from  the  other  two  experimental 
periods  to  validate  the  model.  Additional  simulations  of  KWRF  loading  tests  which  did 
not  reach  steady  state  from  all  three  experimental  periods  were  performed  to  compare  the 
model's  ability  to  predict  clarifier  failure  due  to  solids  overloading  to  that  of  both  limiting 
total  solids  flux  theory  and  the  gravity-flux  constraining  model. 


CHAPTER  2 

ONE-DIMENSIONAL  MODELING  OF  SECONDARY  CLARIFEERS  USING  A 

CONCENTRATION  AND  FEED  VELOCITY-DEPENDENT  DISPERSION 

COEFFICIENT 


Introduction 

The  limiting  total  solids  flux  concept  is  used  for  the  design  of  sludge  thickeners  and 
the  thickening  region  of  activated  sludge  secondary  clarifiers,  where  the  total  solids  flux  is 
the  sum  of  the  solids  flux  due  to  gravity  settling  and  the  solids  flux  due  to  bulk  downward 
movement  of  liquid.  This  concept  originated  with  Coe  and  Clevenger  (1916),  who 
suggested  that  if  a  layer  in  a  suspension  has  a  lower  total  solids-handling  capacity  than  the 
overlying  layer,  it  will  be  unable  to  discharge  solids  as  rapidly  as  they  are  received  and  will 
therefore  grow  in  thickness.  If  a  given  layer  has  a  higher  total  solids-handling  capacity 
than  the  layer  above,  its  thickness  will  decrease  or  remain  infinitesimal.  The  layer  with  the 
lowest  total  solids-handling  capacity  therefore  limits  the  throughput  of  the  thickener.  If 
the  thickener  is  overloaded,  this  layer  (which  contains  the  limiting  solids  concentration) 
will  ultimately  reach  the  liquid  surface.  Yoshioka  et  al.  (1957)  and  Hasset  (1964) 
developed  graphical  procedures  for  computing  the  value  of  the  limiting  total  solids  flux 
that  are  accepted  for  secondary  clarifier  design  (Vesilind,  1968;  WPCF,  1985;  Metcalf  and 
Eddy,  1991).  A  feature  common  to  these  approaches  is  the  postulate  that  the  settling 
velocity  of  sludge  particles  in  the  hindered  settling  regime  is  a  function  only  of  the  local 


suspended  solids  concentration,  as  proposed  by  Kynch  (1952)  in  his  modeling  of  batch 
settling.  He  used  the  method  of  characteristics  to  solve  a  partial  differential  equation 
(PDE)  model. 

A  family  of  one-dimensional,  dynamic  clarifier  models  was  developed  based  on  limiting 
total  solids  flux  theory  (Bryant,  1972;  Tracy  and  Keinath,  1973;  Lessard  and  Beck,  1993). 
These  models  adjust  the  thickness  of  the  layer  with  the  limiting  solids  concentration  so  as 
to  satisfy  the  limiting  total  solids  flux  constraint.  As  a  result,  they  give  steady-state 
concentration  profiles  having  four  distinct  values  in  the  clarifier  (Fig.  2-1).  Above  the  feed 
layer,  the  solids  concentration  is  very  low.  The  feed  layer  has  an  intermediate 
concentration  that  is  less  than  the  feed  concentration.  Below  the  feed  layer,  the  sludge 
blanket  has  a  concentration  equal  to  the  limiting  concentration,  whereas  the  concentration 
of  the  bottom  layer  will  be  higher,  as  set  by  mass  balance.  Petty  (1975)  solved,  using  the 
method  of  characteristics,  a  partial  differential  equation  model  for  the  clarifier  and  raised 
questions  as  to  whether  the  limiting  flux  is  appropriate  in  all  cases. 

A  second  family  of  models  is  based  on  a  modification  of  the  limiting  flux  constraint 
(Stenstrom,  1976;  Hill,  1985;  Vitasovic,  1986,  1989;  Takacs  et  al.,  1991).  Rather  than 
constraining  all  layers  above  the  bottom  layer  to  concentrations  less  than  or  equal  to  the 
limiting  concentration,  they  constrain  only  the  gravity  flux  term.  This  entails  setting  the 
value  of  the  downward  gravity  flux  from  a  given  layer  in  the  thickening  zone  to  the 
minimum  of  the  gravity  flux  calculated  for  that  layer  and  the  gravity  flux  calculated  for  the 
layer  below.  This  approach  avoids  the  necessity  of  calculating  the  limiting  flux  and  limiting 
concentration  and  gives  a  more  realistic  concentration  profile  in  the  thickening  zone. 
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Fig.  2-1.  Example  of  clarifier  concentration  profile  obtained  using  total  limiting  flux 
constraint. 


An  alternative  approach  for  obtaining  a  realistic  concentration  profile  is  to  add  a 
dispersive  flux  term  in  the  mass  balance  for  each  layer  (Anderson  and  Edwards,  1981;  Lev 
et  al.,  1986;  Hamilton  et  al.,  1992).  Adding  a  dispersion  term  converts  the  model  equation 
from  a  hyperbolic  to  a  parabolic  PDE  which  eliminates  problems  with  multiple  solutions 
encountered  using  the  hyperbolic  equation.  Anderson  and  Edwards  (1981)  included  a 
dispersion  term  in  their  model  for  peripheral  feed  clarifiers.  It  is  noteworthy  that  the 
dispersion  coefficient  was  not  constant,  varying  with  position.  Lev  et  al.  (1986)  extended 
the  analysis  of  Petty  (1975)  to  include  the  clarification  zone  and  noted,  as  had  Petty 
(1975),  that  the  limiting  flux  constraint  has  a  limited  range  of  validity  and  that  its 
imposition  could  lead  to  erroneous  results.  They  included  a  dispersion  term  in  a  dynamic 
clarifier  model  and  reported  that  it  yielded  correct  dynamic  behavior.  Hamilton  et  al. 
( 1 992)  modeled  a  pre-denitrification  process,  employing  a  constant  dispersion  coefficient 
in  the  secondary  clarifier  component,  and  proposed  a  method  for  calculating  sludge 
blanket  height. 

Alternative  approaches  emphasize  the  interaction  between  solid  and  liquid  phases  at 
high  solids  concentrations.  Hartel  and  Popel  (1992)  postulated  that  settling  velocity  is 
affected  both  by  underlying  layers  in  the  thickening  zone  and  by  overlying  layers  in  the 
compression  zone,  as  well  as  the  local  suspended  solids  concentration.  They  employed  a 
correction  function  that  reduces  the  settling  velocity  applied  in  the  thickening  zone  based 
on  location  in  the  clarifier  relative  to  the  feed  layer  and  the  position  of  the  compression 
concentration.  They  defined  the  compression  concentration  position  as  the  point  of 
transition  between  hindered  settling  and  compression  and  gave  a  procedure  for  calculating 
it.  In  their  model,  the  gravity  flux  is  the  product  of  the  correction  function,  the  calculated 
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settling  velocity  based  on  concentration,  and  the  layer  concentration.  George  and  Keinath 
(1978)  added  a  liquid  phase  mass  balance  describing  the  change  in  the  upward  velocity  of 
displaced  fluid  with  depth,  and  included  in  their  model  a  settling  velocity  equation  that 
depended  on  the  local  concentration  gradient  as  well  as  the  local  concentration.  It  was 
still  necessary  to  impose  a  limiting  flux  constraint  due  to  the  model's  inability  to  predict 
rising  blankets  under  overloaded  conditions  (Hill,  1985;  George,  1976).  Others  have  also 
questioned  the  Kynch  proposition  that  settling  velocity  depends  solely  on  local  solids 
concentration  (Tiller,  1981;  Fitch,  1983;  Font,  1988).  In  another  approach,  thickening  is 
viewed  in  terms  of  transport  of  mass  and  momentum  in  a  non-rigid  saturated  porous 
medium  (Kos,  1977;  Kos  and  Adrian,  1974;  Landman  et  al.,  1988;  Leonhard,  1993;  Tiller 
and  Hsyung,  1993). 

The  work  of  Takacs  et  al.  (1991),  which  employs  the  gravity  flux  constraint,  is  notable 
in  that  it  presents  an  excellent  match  to  the  full-scale  data  set  of  Pflanz  (1969).  In  the 
present  study,  it  is  shown  that  this  approach  can  be  reinterpreted  as  modeling  with  a 
concentration-dependent  dispersion  coefficient.  The  two  interpretations  give  differing 
results  if  the  number  of  layers  into  which  the  clarifier  is  divided  is  changed.  It  is  shown 
that  for  finer  discretizations,  the  dispersion  interpretation  leads  to  better  fits  with 
experimental  data.  The  expression  for  the  concentration-dependent  dispersion  coefficient 
is  then  simplified  without  decreasing  model  performance.  The  final  model  modification 
was  to  incorporate  a  dependence  of  the  dispersion  coefficient  on  clarifier  feed  velocity. 
The  model  was  applied  to  data  from  experiments  we  conducted  on  a  full-scale  clarifier  at 
the  Kanapaha  Water  Reclamation  Facility  in  Gainesville,  FL. 


Initial  Model  Development 

Takacs  et  al.  (1991)  modified  the  Vesilind  (1968)  equation  for  settling  velocity  to 
account  for  the  fact  that  the  settling  velocity  decreases  as  concentration  approaches  zero. 
They  used  this  equation  in  a  model  employing  the  gravity  flux  constraint.  To  test  their 
model,  they  used  three  cases  of  experimental  data  presented  by  Pflanz  (1969).  Takacs  et 
al.  (1991)  reduced  the  two-dimensional  data  sets  to  one-dimensional  forms  with  respect  to 
depth.  As  each  of  the  reduced  data  sets  involved  ten  depths,  they  modeled  the  clarifier  as 
consisting  often  layers.  Their  model  gave  good  predictions  of  effluent  suspended  solids 
as  well  as  excellent  matches  with  the  solids  concentration  in  the  thickening  zone.  It  was 
found,  however,  that  when  the  number  of  layers  in  their  model  is  increased  to  twenty 
without  changing  model  parameters,  the  model  performance  deteriorates  considerably 
(Fig.  2-2).  Furthermore,  at  finer  discretizations  (e.g.,  50  layers)  model  performance  is 
worse  than  with  a  discretization  of  10  layers,  even  with  parameters  optimally  fitted  for  that 
level  of  discretization  (Fig.  2-3).  Ideally,  model  performance  should  improve  with  an 
increase  in  the  degree  of  discretization.  In  the  following,  the  model  of  Takacs  et  al. 
(1991)  is  analyzed  and  an  alternative  approach  is  developed  that  achieves  this  objective. 

A  mass  balance  on  the  thickening  zone,  using  the  Takacs  et  al.  (1991)  expression  for 
the  gravity  settling  velocity  without  a  flux  constraint,  gives: 

5z  dd  /dt  =  (Qu/A)  Ci.,  -  (Qu  /A)  Ci  +  V^-i  CM  -  Vs>i  C,  (1) 

Vs,,  =  min{V0  [exp(-b(Ci-Cmin))  -  exp(-bp(Ci-Cmin))] ,  V,^  }  (2) 
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Fig.  2-2.  Comparison  of  10-  and  20-layer  versions  of  the  Takacs  et  al.  (1991)  model  using 
parameters  estimated  for  the  10-layer  version  of  the  model  applied  to  case  1  of  the 
Pflanz  data. 
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Fig.  2-3.  Comparison  of  10-and  50-layer  versions  of  the  Takacs  et  al.  (1991)  model,  with 
optimal  parameters  estimated  for  each  version  applied  to  case  1  of  the  Pflanz  data. 
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where  C;  is  the  \th  (from  the  top)  layer  suspended  solids  concentration,  Cmu,  is  the 
nonsettleable  suspended  solids  concentration,  V0  and  b  are  the  Vesilind  (1968)  settling 
parameters,  bp  is  a  settling  parameter  characteristic  of  low  suspended  solids 
concentrations,  Vmax  is  the  highest  settling  velocity  achieved  by  sludge  floes,  Vs,;  is  the 
gravity  settling  velocity  from  layer  i,  Qu  is  the  underflow  flow  rate,  A  is  the  clarifier  cross- 
sectional  area,  and  8z  is  the  layer  thickness.  Incorporating  the  constraint  on  gravity  flux 
(i.e.,  setting  the  value  of  the  downward  gravity  flux  from  a  given  layer  to  the  minimum  of 
the  gravity  flux  calculated  for  that  layer  and  the  gravity  flux  calculated  for  the  layer 
immediately  below)  modifies  eq.  (1)  to: 

5z  dCi  /dt  =  (Qu  /A)  Ci.,  -  (Qu  /A)  Ci  +  min[Vs,i.,  CM,  VSii  CJ  -  min[VSii  C;,  Vs,i+i  Cm]     (3) 

An  examination  of  the  Takacs  model  as  applied  to  the  Pflanz  data  shows  that  the 
constraint  on  the  gravity  flux  becomes  active  at  a  certain  layer  in  the  thickening  zone  and 
remains  active  throughout  all  lower  (i.e.,  with  higher  i)  layers.  In  that  region,  the  material 
balance  for  each  layer  becomes: 

6z  dCi  /dt  =  (Qu  /A)  Cm  -  (Q„  /A)  Q  +  VtfQ  -  Vs,i+1Ci+1  (4) 

An  alternate  approach  is  to  add  a  dispersion  term  to  eq.  (1),  in  which  case  the  mass 
balance  for  layer  i  is: 

5z  dC,  /dt  =  (Qu  /A)  Cm  -  (Qu  /A)  C;  +  V^iQ.,  -  Vs,iCi  -  Dfc.y  (Ci-Ci.,)/5z 

+  Di,i+1(Ci+1-C,)/8z  (5) 
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where  Dm,;  is  the  dispersion  coefficient  for  the  dispersive  flux  from  layer  i  to  layer  i-1  and 
Di,i+i  the  coefficient  for  the  flux  from  layer  i+1  to  layer  i.  Eq.  (5)  becomes  identical  to  eq. 
(4)  if  one  uses  the  following  concentration-dependent  expression  for  the  dispersion 
coefficient: 

Di,i+1(Ci,Ci+1)  =  8z  {Vs,iCi  -  Vs,i+1Ci+1  }/(Ci+,  -  CO  (6) 

The  agreement  with  the  gravity  flux  constraint  is  complete  if  one  imposes  the  physical 
constraint 

rWCfcCw)  >  o  (7) 

This  is  because  the  expression  for  the  dispersion  coefficient  becomes  negative  exactly 
when  the  gravity  flux  constraint  becomes  inactive.  With  the  above  constraint  applied  to 
the  dispersion  coefficient,  eq.  (5)  becomes  equivalent  to  eq.  (3)  throughout  the  clarifier. 

Eq.  (6)  implies  that  the  dispersion  coefficient  disappears  as  the  layer  thickness  5z  — >  0. 
This  is  not  physically  correct.  In  terms  of  the  original  limiting  gravity  flux  formulation,  as 
5z  -»  0,  eq.  (3)  takes  the  form  of  eq.(l)  throughout  the  clarifier,  i.e.,  the  constraint  on  the 
gravity  flux  disappears  everywhere. 

To  correct  this  problem,  eqs.  (6)  and  (7)  are  modified  to 

Du+,(Ci,Ci+1)  =  max[a(Vi.iCi  -  Vs,i+iCi+1)]/(Ci+1  -  Q)  ,  0]  (8) 

i.e.,  the  new  parameter  a  replaces  5z.  In  this  case  as  5z  — >  0,  the  dispersion  coefficient  D 
converges  to  the  finite  quantity 
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D(C)=max{-a^^,0}  (9) 

dC 


If  the  number  of  layers  into  which  the  clarifier  is  subdivided  changes,  so  too  will  the 
layer  thickness  5z.    In  the  model  with  dispersion  (eq.  (8)),  a  will  remain  constant, 
whereas  the  constrained  gravity  flux  approach  (eq.  (3))  is  equivalent  to  modifying  a 
according  to  the  change  in  8z.  Thus  the  two  approaches  will  provide  different  results  if 
the  level  of  discretization  is  changed.  Figure  2-4  presents  concentration  profiles  obtained 
by  each  approach  for  the  Pflanz  data  (first  case)  with  a  discretization  of  20  layers  using  the 
parameters  that  Takacs  et  al.  (1991)  reported  for  the  10-layer  discretization.  It  is  clear 
that  eq.  (8)  provides  a  better  fit.  Cases  2  and  3  of  the  Pflanz  data  provide  similar  results 
(data  not  shown).  It  is  concluded  from  this  result  that  the  Takacs  et  al.  model  would 
benefit  by  removing  its  dependence  on  the  level  of  discretization.  This  uncoupling  is 
accomplished  by  recasting  the  equations  resulting  from  the  gravity  flux  constraint  as  a 
dispersive  flux  term,  as  presented  in  eq.  (8). 

It  would  be  advantageous  to  simplify  the  expression  for  the  concentration-dependent 
dispersion  coefficient  (eq.  (8)).  Since  Cmin,  exp(-bp  (Ci-Cmin),  and  the  Vmax  constraint  are 
only  significant  in  layers  of  low  concentration  where  the  first  argument  of  the  max  operator 
in  eq.  (8)  is  negative  and  therefore  not  used,  these  terms  can  be  neglected,  yielding 

Du+i(Ci,Ci+1)  =  max{oc  V„  [exp(-b  Q)  Q  -  exp(-b  Ci+1)  Ci+1]/(Ci+1  -  Q),  0}  (10) 

If  Ci+i  is  close  to  C; ,  exp(-b  C;+i)  C;+i  and  exp(-b  C;)  Cj  can  be  approximated  by  their  first 
order  Taylor  series  expansions  about  the  geometric  mean  of  the  concentrations 
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Fig.  2-4.  Comparison  of  concentration  profiles  obtained  by  20-layer  versions  of  the 
Takacs  et  al.  (1991)  model  (eq.  (3))  and  the  model  with  dispersion  (eq.  (8))  to 
experimental  data  (case  1  of  the  Pflanz  data).  Parameters  reported  for  the  10-layer 
model  of  Takacs  et  al.  (1991)  were  employed  to  generate  model  fits. 
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( Cii+1=A/Ci  Ci+1 ).  In  this  case,  eq.  (10)  becomes 

Di,i+1(Ci,C1+1)  -  max{a  V0  (bCu+1-l)  exp(-b  Cu+1),  0}  (11) 

The  first  argument  of  the  max  operator  in  the  above  expression  becomes  negative  for 
Cii+1  <  1/b  and  attains  a  maximum  when  Ci?j+1=  2/b.  It  makes  physical  sense  that  at  high 

concentrations  (>2/b)  the  dispersion  coefficient  decreases  with  increasing  concentration, 
as  eq.  (1 1)  implies.  It  does  not,  however,  make  physical  sense  that  the  dispersion 
coefficient  decreases  as  the  concentration  decreases  at  low  concentrations  (<  21b),  as  eq. 
(1 1)  also  implies.  Therefore  Dij+i(Ci,C;+i)  is  set  equal  to  its  maximum  (=  a  V0  exp(-2)) 
for  concentrations  C;  i+1  <  2/b.  The  expression  for  the  dispersion  coefficient  now 
becomes: 


Du+,(C„C1+1)  = 


a  V0  (bCu+1  - 1)  exp(-b  Cu+1 )  for  Cu+1  >  2  /  b 
aV0exp(-2)  forC11+,  <2/b 


The  above  expression  involves  parameters  from  the  settling  equation  (Vo,  b)  in  addition  to 
the  introduced  parameter  a.  The  dispersion  parameters  can  be  decoupled  from  the  settling 
parameters  by  introducing  P  =  b,  Ccrit  =  2/b,  and  Dmax  =  a  V0  exp(-2),  in  which  case  eq. 
(12)  is  replaced  by: 


JDmax[l  +  3(Cu+1  -Ccm)]exp[-P(Cu+1  -Ccnt)]  for  CM+1  >  C^ 
I  Dmax  for  Cu+1  <  Ccnt 
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Fitting  P  and  Ccrit  instead  of  computing  them  from  b  gives  more  degrees  of  freedom  and 
therefore  potentially  better  fits,  but  at  the  expense  of  having  extra  estimated  parameters. 

The  above  expression  (eq.  (13))  for  the  dispersion  coefficient  D  is  the  simplest 
function  that  has  the  following  features: 

•  Sets  D  equal  to  a  constant  Dmax  for  low  concentrations  (less  than  Ccrit) 

•  Decreases  D  exponentially  with  increasing  C  for  high  concentrations.  The  physical 
justification  for  this  is  that  viscosity  increases  with  increasing  suspended  solids 
concentration. 

•  Provides  for  a  smooth  transition  between  the  constant  D  region  and  the  exponential 
decay  region.  The  contribution  of  the  factor  l+P(C-Ccrit)  is  to  eliminate  a 
discontinuity  (corner)  in  the  slope  at  C  =  Ccrit . 

The  complete  model  equations  as  applied  to  a  cylindrical  clarifier 
(Fig.  2-5)  are  now  presented,  with  Qf  denoting  the  flow  rate  into  the  clarifier,  Qe  the 
effluent  flow  rate,  Qu  the  underflow  (return  activated  sludge)  flow  rate,  A  the  cross- 
sectional  area,  Uqthe  overflow  velocity  (Qe/A),  and  Ub  the  underflow  velocity  (Qu/A) . 
The  clarifier  is  subdivided  into  layers  of  thickness  8z,  with  numbering  from  top  to  bottom. 
In  this  geometry,  the  effluent  concentration  (Ce)  will  be  the  concentration  of  the  first 
layer,  whereas  the  return  activated  sludge  (RAS)  concentration  (Cu)  will  be  the 
concentration  of  the  bottom  layer. 

For  the  top  layer  (i  =  1): 

5z  dCi/dt  =  Uq  C2  -  UqCi  -  VMCi  +  D,,2  (C2  -  Ci)/  8z  (14) 


Of 
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Fig.  2-5.  Clarifier  geometry  used  for  initial  model  development.  Clarifier  is  divided  into  n 
layers. 
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For  the  \th  layer  in  the  clarification  section: 

5z  dCi/dt  =  Uq  Ci+,  -  UqCi  +  Vs>i.,Ci.i  -  Vs,;Ci  -  D;.,,j  (Cj-CM)/8z  +  D;,i+1  (Ci+1-Ci)/8z     (15) 

For  the  layer  receiving  clarifier  feed: 

5z  dCi/dt  =  (Qf /A)  Cf  -  Uq  Cj  -  UbCi  +  V^Q.,  -  Vs,iCi  -  D^u  (Q-Cj-O/Sz 

+  Du+1  (Ci+i-Ci)/5z  (16) 

For  the  \th  layer  in  the  thickening  zone: 

5z  dCi/dt  =  Ub  Cm  -  Ub  Q  +  V^iCu  -  Vs,iCi -  Dj.jj  {CrCx.x)lbz  +  Du+1  (Ci+i-Ci)/5z    (17) 

And  finally  for  the  bottom  layer: 

8z  dd/dt  =  Ub  CM  -  UbCi  +  Vs.i.jCi.i  -  D{.hi  (Q-CwySz  (18) 

The  above  equations  are  essentially  the  result  of  discretizing  the  parabolic  partial 
differential  equation 


f-f-f^of' 


where  U  =  Uq  in  the  clarification  section,  U  =  -  Ub  in  the  thickening  zone,  and  Gs  is  the 
gravity  settling  flux. 
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Application  of  Clarifier  Model  to  the  Pflanz  Data 

Table  2-1  presents  the  fit  of  the  previously  presented  model  to  the  three  sets  of  Pflanz 
data  as  modified  by  Takacs  et  al.  (1991).  Parameters  were  estimated  using  the  Levenberg- 
Marquardt  algorithm  (Marquardt,  1963;  Press  et  al.,  1989;  Cuthbert,  1987).  The 
objective  function  for  parameter  estimation  was  the  sum  of  the  squares  of  the  relative 
errors  (SSRE)  between  observed  and  model  concentrations.  As  was  done  by  Takacs  et  al. 
(1991),  a  separate  set  of  parameter  values  was  estimated  for  each  case. 

The  clarifier  was  discretized  into  50  layers  and  the  feed  layer  was  set  in  a  position 
consistent  with  that  chosen  by  Takacs  et  al.  (1991)  in  their  10-layer  discretization.  The 
table  shows  the  fit  with  the  dispersion  expression  of  eq.  (13)  for  three  cases: 

•  One  dispersion  parameter  estimated  (Dmax).  The  other  two  parameters  are  calculated 
from  settling  parameters  as  implied  by  eq.  (12),  i.e.,  P  =  b  and  Ccrit  =  2/b.  We  refer  to 
this  as  the  Dmax  model. 

•  Two  dispersion  parameters  estimated  (Dmax  and  Ccrjt)  with  P  =  b.  We  refer  to  this  as 
the  Dmax-Ccnt  model. 

•  Three  dispersion  parameters  estimated  (Dmax,  Ccrit  and  P).  We  refer  to  this  as  the  D,™- 
Ccrit-P  model. 

As  can  be  seen  from  Table  2-1,  the  D^  model  did  quite  well  (SSREs  0.087,  0.067 
and  0.077  for  the  three  respective  cases  of  the  Pflanz  data)  in  relation  to  the  10-layer 
Takacs  et  al.  (1991)  model  (SSREs  0.275,  0.254  and  0.157).  Recall  that  the  Takacs  et  al. 
(1991)  model  performed  worse  for  a  50-layer  discretization  (Fig.  2-3).  Some 
improvement  is  obtained  by  the  Dmax-Ccrit  model  (SSREs  0.085,  0.060  and  0.056),  while 
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the  Dmax-CcrirP  model  offers  no  further  improvement.    The  estimated  parameters  are 
given  in  Table  2-2.  A  10-layer  discretizaton  of  the  model  with  dispersion  (eq.  (13))  was 
also  investigated  and  gave  SSREs  of  0.258,  0.090,  0.187  when  the  Dmax  model  was  used 
and  0.250,  0.089,  0. 185  when  the  Dmax-Ccrit  model  was  employed  (data  not  shown).  As 
expected,  these  SSREs  were  higher  than  obtained  with  the  50-layer  discretization. 

Experimental  Measurements  on  a  Full-Scale  Clarifier 

Field  tests  were  conducted  at  the  Kanapaha  Water  Reclamation  Facility  (KWRF)  to 
collect  data  for  testing  the  clarifier  model.  The  KWRF  is  a  38,000  m3/d  (10  Mgal/d)  plant 
utilizing  the  Ludzack-Ettinger  process  for  nitrification,  denitrification,  and  carbon 
oxidation  (Ludzack  and  Ettinger,  1962).  Clarifier  loading  tests  were  conducted  using  one 
of  the  plant's  four  secondary  clarifiers.  Batch  sludge  settling  tests  were  carried  out  in 
parallel  with  the  loading  tests  to  provide  data  for  determining  settling  equation  parameters. 

Clarifier  Test  Procedure 

The  secondary  clarifiers  were  28.96  m  (95  ft)  in  diameter  with  a  3.66  m  (12  ft) 
sidewall  depth  and  a  4.6  m  (15  ft)  depth  at  the  center.  The  influent  entered  at  a  central 
feedwell  that  was  bounded  by  an  annular  baffle  (shroud)  extending  from  above  the  liquid 
surface  to  2.44  m  (8  ft)  below  the  liquid  surface.  Effluent  left  the  system  over  peripheral 
and  radial  weirs.  The  clarifier  feedwell  occupied  approximately  28  percent  of  the  cross- 
sectional  area.  RAS  was  removed  continuously  via  a  rotating  multiple-pipe  suction  system 
with  four  draw-offs  approximately  0.23  m  (9  inches)  above  the  bottom  and  spaced  at 
distances  of  1.5,  5.2,  8.6,  and  12.2  m  from  the  center  of  the  clarifier.  Waste  activated 
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sludge  (WAS)  was  removed  by  periodic  pumping  from  a  central  sump.  During  the  loading 
tests,  the  influent  flow  rate  to  the  clarifier  was  controlled  by  flooding  the  distribution  box 
that  fed  the  four  clarifiers  and  then  adjusting  the  flow  rate  manually  using  the  in-line  valves 
between  the  distribution  box  and  the  four  clarifiers.  Flow  adjustments  were  made  with 
reference  to  a  hand-held  meter  that  displayed  the  effluent  flow  rate  of  the  test  clarifier. 
The  underflow  rate  was  controlled  by  adjusting  the  speed  of  the  RAS  pump  with  reference 
to  a  flow  meter  on  the  RAS  line. 

A  loading  test  was  initiated  by  setting  the  influent  and  underflow  flow  rates  to  selected 
values  and  typically  lasted  8-10  hours.  The  sludge  blanket  height  was  measured  at  15 
minute  intervals  using  a  5  cm  ID  transparent  plastic  tube  (Sludge  Judge).  In  loading  tests 
where  a  steady-state  was  achieved  as  judged  from  the  blanket  height,  data  collection  was 
continued  for  at  least  two  hours  into  the  steady- state  period.  In  tests  where  a 
continuously  rising  blanket  was  observed,  data  collection  was  ended  when  the  sludge 
blanket  approached  the  effluent  weirs.  Mixed  liquor  flowing  into  the  clarifier,  clarifier 
effluent,  RAS  and  WAS  were  sampled  hourly:  mixed  liquor  starting  at  the  beginning  of 
the  test  and  the  remaining  streams  beginning  at  the  time  that  the  sludge  blanket  reached  a 
stable  height  or  when  it  became  apparent  that  the  sludge  blanket  would  not  stop  rising. 
Samples  were  stored  on  ice  until  analysis  for  suspended  solids,  which  was  performed  on 
the  same  day  as  the  loading  test.  Concentration  profiles  were  determined  in  selected  tests 
by  collecting  samples  at  0.61  m  intervals  through  the  clarifier  depth.  The  sampling 
apparatus  allowed  collection  of  samples  at  three  depths  simultaneously.  Concentration 
profiles  were  measured  at  one  point  inside  the  shroud  and  another  point  halfway  between 
the  shroud  and  the  peripheral  wall  of  the  clarifier. 
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Batch  Settling  Tests 

A  water-bath  enclosed  six-column  settling  apparatus  was  constructed  for  the  batch 
sludge  settling  tests  based  on  the  design  of  Wahlberg  and  Keinath  (1986).  Tests  were 
carried  out  at  sludge  concentrations  ranging  from  2000  to  14000  g/m  .  The 
concentrations  were  achieved  by  mixing  RAS,  mixed  liquor,  and  clarifier  effluent  in 
selected  ratios.  Samples  of  the  secondary  effluent,  mixed  liquor,  and  RAS  collected  for 
the  settling  test  were  taken  for  later  suspended  solids  analysis.  Sludge  in  the  columns  was 
mixed  for  5  minutes  immediately  prior  to  the  settling  test  using  compressed  air  introduced 
through  air  stones.  After  mixing,  the  interface  height  was  recorded  every  two  minutes 
until  the  compression  phase  was  reached.  Total  test  duration  ranged  from  30  minutes  to 
2.5  hours  depending  on  sludge  concentration. 

The  settling  velocity  (Vs)  at  each  initial  suspended  solids  concentration  (C)  was 
determined  from  the  slope  of  data  points  lying  along  the  initial  linear  portion  of  the 
interface  height  versus  time  curve.  The  expression  (V,  =  Vo  e"*0  )  (Vesilind  1968)  was 
used  in  finding  the  Vesilind  parameters  V0  and  b  by  least  squares  linear  regression  on  the 
logarithms  of  the  settling  velocities  and  the  corresponding  sludge  concentrations,  as 
recommended  by  Daigger  (1995).  (A  settling  test  refers  to  the  six  settling  trials  carried 
out  simultaneously  in  the  multi-column  apparatus.)  A  total  of  seven  settling  tests  were 
carried  out  over  the  three  week  experimental  period.  Values  of  V0  and  b  were  found  for 
each  test  and  averaged.  One  settling  test  gave  values  of  V0  and  b  that  were  more  than  two 
standard  deviations  from  the  mean  of  all  seven  V0  and  b  values.  Data  from  this  test  were 
rejected,  then  velocity  versus  concentration  data  from  the  six  accepted  tests  were  pooled 
(Fig.  2-6)  and  a  linear  regression  was  performed  to  obtain  a  single  set  of  Vesilind 
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Fig.  2-6.  Pooled  data  from  batch  settling  tests  after  discarding  outliers.  Model  line 
represents  fit  of  Vesilind  equation  with  V0  =  182.9  m/d  and  b  =  0.3055  m3/kg. 
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parameters  for  the  test  period  (V0  ■  182.9m/d  and  b  =  0.3055  m3/kg). 


Additional  Measurements 

The  60-minute  nonsettleable  suspended  solids  concentration  (Method  2540F,  APHA 
et  al.,  1992)  was  measured  in  duplicate.  A  sample  of  mixed  liquor  was  settled  for  60 
minutes  in  a  plastic  one  L  graduated  cylinder,  then  a  supernatant  volume  of  150  mL  was 
withdrawn  for  subsequent  TSS  analysis.  The  mean  of  the  60-minute  nonsettleable 
suspended  solids  concentrations  was  3.2  g/m3.  This  was  higher  than  the  TSS  of  a 
significant  number  of  clarifier  effluent  grab  samples.  The  mean  of  these  samples  (2.4 
g/m3)  was  therefore  used  as  the  nonsettleable  suspended  solids  concentration  for 
simulations.  Samples  for  total  suspended  solids  analysis  were  filtered  through  glass  fibre 
filters  having  an  average  pore  size  of  1 .2  jim  (Whatman  GF/C).  Filtered  residues  were 
dried  to  constant  weight. 

Application  of  Clarifier  Model  to  KWRF  Data  and  Further  Model  Development 

The  clarifier  model  presented  in  the  previous  section  assumes  that  the  clarifier  can  be 
regarded  as  a  cylinder  with  underflow  removed  uniformly  from  the  bottom.  In  the  KWRF 
clarifier,  underflow  is  removed  by  a  hydraulic  suction  system  which  has  intake  pipes 
spaced  along  the  sloping  floor  of  the  bottom,  conical  section.  Because  the  pipe  intakes  are 
at  different  depths,  the  concentrations  of  sludge  withdrawn  at  the  different  depths  will  be 
substantially  different,  and  therefore  RAS  removal  cannot  be  modeled  as  if  it  were 
withdrawn  from  a  single  layer.    The  clarifier  model  was  therefore  modified  by  including  a 
conical  section  at  the  bottom  of  the  main  cylindrical  section  (Fig.  2-7)  The  conical  section 
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Fig.  2-7.  Schematic  diagram  of  full-scale  secondary  clarifier  at  the  Kanapaha  Water 
Reclamation  Facility,  showing  shroud,  bottom  conical  section  and  model 
discretization. 
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is  divided  into  p  layers,  whereas  the  cylindrical  section  is  divided  into  n  layers.  In  each  of 
the  layers  of  the  conical  section,  sludge  is  withdrawn  at  a  rate  equal  to  the  change  in 
cross-sectional  area  from  the  top  to  the  bottom  of  the  layer  multiplied  by  the  underflow 
velocity: 

Qr,i  =  (A-A+i)Ub  (20) 

This  results  in  a  constant  underflow  velocity  throughout  the  conical  section.  The  total 
RAS  flow  rate  (Q„)  is  the  summation  of  the  withdrawals  from  the  layers  of  the  conical 
section: 

n+p 

Qu=  I  Q  (21) 

The  RAS  concentration  was  calculated  as  the  flow- weighted  average  of  the  concentrations 
of  layers  in  the  conical  section. 

The  model  was  also  modified  to  account  for  the  presence  of  a  shroud  in  the  upper 
section  of  the  clarifier.  The  cross-sectional  area  of  the  clarifier  available  for  overflow  in 
the  upper  region  of  the  clarifier  (i.e.,  from  0  to  2.44  m  below  the  water  surface)  is  72%  of 
the  cross-section  below  the  shroud  (Fig.  2-7). 

Based  on  observations  of  density  current  flows  in  prototype  scale  clarifiers  (Andersen, 
1945),  the  feed  in  the  clarifier  model  was  input  to  the  layer  above  the  first  layer  having  a 
concentration  greater  than  the  feed  concentration.  This  was  accomplished  by  the 
following  recursive  procedure.  A  position  for  the  feed  layer  was  initially  assumed.  The 
concentration  profile  was  then  calculated  and  used  to  update  feed  layer  position.  The 
latter  two  steps  were  repeated  until  convergence  was  achieved.  In  rare  cases  the  above 
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procedure  did  not  converge  to  a  single  layer  and  instead  began  to  oscillate  between  two 
adjacent  layers.  In  those  cases,  the  higher  of  the  two  layers  was  chosen  as  the  feed  layer. 

Data  from  nine  clarifier  loading  tests  at  the  KWRF  in  which  a  steady-state  blanket  was 
attained  and  the  concentration  profile  was  measured  were  used  to  calibrate  and  evaluate 
the  model.  These  tests  involved  a  range  of  underflow  and  overflow  rates  (Table  2-3). 
During  the  same  experimental  period,  five  other  tests  failed  to  reach  a  steady  state  because 
of  overloading  (i.e.,  continually  rising  blanket)  and  in  one  other  test  a  steady  state  was 
reached  but  the  concentration  profile  was  not  measured.  An  example  of  a  loading  test  in 
which  a  steady  blanket  level  was  achieved  is  shown  in  Figure  2-8a,  whereas  Figure  2-8b 
shows  an  example  of  a  test  in  which  the  blanket  continued  rising  throughout  the 
experimental  period. 

Model  parameters  were  determined  as  follows:  V0,  b,  and  Cmi„  were  obtained  by 
analysis  of  batch  settling  data  and  clarifier  effluent  samples,  as  described  previously.  The 
parameters  bp,  Vmax,  and  Dmax  were  obtained  by  least-squares  nonlinear  regression  on 
concentration  profiles  as  explained  below.  The  remaining  parameters  (Ccrit  and  P)  were 
either  estimated  by  nonlinear  regression  on  concentration  profiles  or  calculated  from  the 
experimentally  determined  b  (P  =  b  and  Ccrit  =  2/b).  The  effluent  (overflow)  and  RAS 
concentrations  used  in  fitting  were  the  average  values  over  the  period  when  the  system 
was  determined  to  be  at  steady-state,  whereas  the  remaining  points  in  the  concentration 
profile  were  from  a  single  set  of  measurements  taken  near  the  end  of  each  loading  test. 

To  determine  model  parameters  for  the  Pflanz  data  (Takacs  et  al.,  1991),  the  sum  of 
the  squares  of  relative  errors  in  concentration  was  used  as  the  performance  measure.  This 
was  selected  because  Takacs  et  al.  (1991)  reported  the  quality  of  their  fits  in  terms  of 
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Fig.  2-8.  Results  of  typical  clarifier  loading  tests:  (a)  test  in  which  a  steady  blanket  level 
was  achieved  (case  4),  (b)  test  in  which  blanket  continued  rising  throughout 
experimental  period  (case  13). 
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relative  errors  in  concentration.  An  alternative  performance  measure  is  the  sum  of  the 

squares  of  the  errors  in  the  logarithms  of  concentrations  (SSELC).  This  measure  provides 

for  better  fits  in  the  thickening  (high  concentration)  section  at  the  expense  of  somewhat 

worse  fits  in  the  clarifying  section  (low  concentration).  As  the  relative  errors  in 

concentration  measurements  are  higher  at  low  concentrations,  this  is  a  desired  tradeoff. 

Therefore  parameter  fitting  with  the  KWRF  data  was  based  on  the  SSELC. 

The  dashed  lines  in  Figure  2-9  show  the  Dmax  model  fit  to  the  KWRF  data.  These 

curves  were  generated  by  a  single  set  of  model  parameters  for  all  nine  cases.  The  fits  in 

some  of  the  cases  (cases  1,  2,  4,  6)  approximated  quite  well  smooth  curves  that  could  be 

drawn  to  represent  the  measured  data.  In  other  cases,  particularly  at  high  loadings  (cases 

7,  8  and  9),  the  fits  are  poor.  The  total  SSELC  for  the  Dmax  model  (Table  2-4)  was  29.2. 

The  Dmax-CCrit  model  and  the  Dmax-CcrirP  model  gave  somewhat  better  fits  overall 

(SSELCs  of  24.9  and  21.8,  respectively),  but  still  performed  poorly  in  cases  of  high 

loading  (Table  2-4).  If,  however,  D^x  is  allowed  to  vary  from  case  to  case,  satisfactory 

fits  can  be  obtained  for  all  cases,  even  with  the  D^x  model  (Table  2-4,  last  column).  This 

leads  to  the  investigatation  of  a  possible  dependence  of  Dmax  on  clarifier  loading.  There  is 

physical  justification  for  correlating  dispersion  to  velocity.  For  example,  Taylor  in  his 

classic  paper  (Taylor,  1953)  reports  a  quadratic  dependence.  Figure  2-10  shows  the 

individually  fitted  Dmax  versus  feed  velocity  Vf  (=  Qf /A).  It  is  observed  that  the  equation 

f      D.+ytVf-Vf,!)2  ifVf>Vf,i 
Dmax=    \  (22) 

I      D,  if  Vf<  Vu 
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Table  2-4.  Results  for  models  with  Dmax  constant  across  all  cases  and  for  the  D,,^  model  with 
Dmax  fitted  individually  for  each  case 


Sum  of  Squares  of  Errors 

in  Logarithms  of  Concentrations 

Case 

D,™  model 

Dnux-Ccht  model 

L'max 

-CcrirP  model 

Dmax  model  with 

individually  fitted  Dmax 

1 

0.882 

1.618 

2.314 

1.084 

2 

0.323 

0.803 

1.610 

0.148 

3 

4.347 

1.783 

0.812 

0.126 

4 

0.245 

0.331 

1.212 

0.073 

5 

1.254 

0.799 

0.440 

0.327 

6 

0.757 

0.547 

0.525 

0.148 

7 

4.381 

3.577 

1.921 

0.385 

8 

1.450 

0.560 

0.184 

0.110 

9 

15.533 

14.918 

12.754 

1.274 

Total 

29.172 

24.935 

21.772 

3.675 
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fits  the  relationship  between  Dmax  and  Vf  quite  well.  The  model  sum  of  errors  squared  is 
only  6%  of  the  total  variation  in  the  fitted  D^  values. 

Two  versions  of  our  clarifier  model  which  incorporated  the  feed-velocity  dependence 
were  further  investigated:  one  in  which  P  and  Ccrit  are  calculated  from  settling  parameters 
and  one  in  which  (3  and  CCTn  are  found  by  model  fitting.  These  are  referred  to  as  the  feed- 
velocity-dependent  Dmax  (FVDDmax)  and  FVDDmax-CcrirP  models,  respectively.  Fits  of 
the  FVDDmax  and  FVDDmax-Ccrit-P  models  to  the  KWRF  data  were  carried  out  using 
SSELC  as  the  objective  function  (see  Table  2-5  for  parameter  values).  These  fits  are 
shown  by  the  dotted  (FVDDmax)  and  solid  (FVDDmax-Ccrit-p)  lines  in  Figure  2-9.  The  fits 
with  these  models  are  greatly  improved  in  most  of  the  cases  over  those  with  constant  Dmax 

In  most  cases,  the  FVDDmax  model  gave  results  almost  indistinguishable  from  those 
obtained  with  the  FVDDmax-Ccrit-P  model.  Table  2-6  compares  profiles  from  the  FVDDmax 
model  to  measured  data  and  gives  values  for  measured  and  predicted  RAS  concentrations. 

RAS  concentrations  were  less  than  concentrations  at  the  bottom  of  the  clarifier  (i.e.,  layer 
50)  because  RAS  was  drawn  off  at  depths  throughout  the  bottom  conical  section.  The 
total  SSELC  for  this  fit  was  5.05.  The  SSELC  for  the  fit  of  the  FVLMWCcrit-p  model 
was  4.34  (data  not  shown).  For  comparative  purposes,  a  fit  was  also  carried  out  on  the 
KWRF  steady-state  data  using  the  Takacs  et  al.  (1991)  model,  with  the  feed  layer 
determined  according  to  the  recursive  procedure  outlined  earlier,  a  conical  section  at  the 
bottom,  and  50  layers.  This  model  gave  a  SSELC  of  45.8. 

The  predictive  capability  of  the  FVDDmax  and  FVDDmax-Ccrit-P  models  was  evaluated 
by  comparing  their  output  to  the  results  of  the  extra  six  cases  from  the  experimental 
period,  none  of  which  was  used  in  model  development.  Both  the  FVDDmax  and 
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Table  2-5.  Parameters  resulting  from  FVDDmax  and  FVDDmax-Ccnt-3  fit  across  nine 
KWRF  cases  for  which  steady  blanket  levels  were  obtained 


Parameters 

Model 

FVDD,^ 

FVDDmax-Ccnl-b 

'max 

173.5 

171.9 

m/d 

bP 

28.8 

27.2 

m3/kg 

D, 

2.95 

3.477 

m2/d 

8 

0.0446 

0.0507 

d 

v«, 

34.7 

37.0 

m/d 

^crit 

* 

11654 

g/m3 

b 

* 

0.00287 

m3/kg 

*  calculated  from  parameter  b 
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FVDDmax-Ccrit-P  models  correctly  predicted  failure  for  the  five  cases  in  which  a  continually 
rising  blanket  was  observed  experimentally  and  also  predicted  success  for  the  "tenth" 
steady-state  case  for  which  a  concentration  profile  was  not  available.  Table  2-7  compares 
the  predictions  of  the  FVDDmax  and  FVDDmax-CcrjrP  models  for  all  15  cases  to  the 
experimental  observations  and  also  to  predictions  of  the  Takacs  et  al.  (1991)  and  total 
limiting  solids  flux  (e.g.,  Coe  and  Clevenger,  1916)  models.  Notably,  the  latter  two 
models  predicted  success  in  all  of  the  cases  where  failure  was  observed.  White  (1976)  and 
Ekama  and  Marais  (1986)  observed  that  total  limiting  flux  theory  can  overpredict  clarifier 
thickening  capacity  by  as  much  as  20%.  Reduction  of  the  calculated  limiting  flux  by  20% 
would  classify  cases  1 1  and  12  as  overloaded,  case  13  as  borderline,  and  cases  14  and  15 
as  still  underloaded.  In  the  latter  two  trials,  the  underflow  velocity  (Ub)  exceeded  the 
critical  underflow  velocity,  i.e.  the  plot  of  total  solids  flux  vs.  concentration  had  no 
minimum,  only  an  inflection  point.  In  these  cases,  the  limiting  total  solids  flux  was 
computed  according  to  the  procedure  used  by  White  (1976)  for  such  situations. 

Conclusions 

As  a  result  of  the  modeling  and  experimental  work  carried  out  in  this  study,  the 
following  conclusions  can  be  drawn: 
•     Although  clarifier  models  incorporating  a  constraint  on  gravity  flux  can  provide 

excellent  simulation  of  experimental  concentration  profiles,  the  flux  constraint 

effectively  disappears  as  the  level  of  discretization  is  increased 
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•  The  gravity  flux  constraint  can  be  recast  as  a  concentration-dependent  dispersion  term, 
which  improves  the  ability  of  the  model  to  fit  experimental  data  as  the  level  of 
discretization  increases 

•  A  further  advantage  of  the  dispersion  model  is  that  the  effect  of  feed  velocity  on 
clarifier  thickening  performance  can  be  accounted  for 

•  Total  limiting  solids  flux  theory  as  well  as  models  incorporating  a  gravity  flux 
constraint  can  fail  to  predict  overloading  of  clarifiers. 

•  To  be  most  valuable,  the  clarifier  model  should  be  integrated  with  a  model  of  the 
activated  sludge  process  under  investigation,  using  for  example  the  IAWPRC 
Activated  Sludge  Model  No.  1  (Henze  et  al.,  1986).  In  this  way,  the  integrated  model 
can  be  employed  to  simulate  the  impacts  of  varying  flow  rates  and  feed  compositions 
on  both  biochemical  and  sludge  thickening  performance. 


CHAPTER  3 
CALIBRATION  OF  A  ONE-DIMENSIONAL  CLARIFIER  MODEL  USING  SLUDGE 

BLANKET  HEIGHTS 


Introduction 

There  has  been  considerable  work  on  developing  easy-to-use  models  for  secondary 
clarifiers.  One  group  of  models  is  based  on  application  of  limiting  solids  flux  theory 
(Tracy  and  Keinath,  1973;  Bryant,  1972;  Lessard  and  Beck,  1993).  Another  approach  to 
clarifier  modeling  is  to  impose  a  constraint  not  on  the  total  solids  flux  but  only  on  the 
gravity  settling  flux  term  (Stenstrom,  1976;  Vitasovic,  1989;  Takacs  et  al.,  1991).  This 
approach  can  yield  a  realistic  solids  concentration  profile  with  respect  to  depth.  An 
alternative  approach  to  the  above  models  is  the  introduction  of  a  dispersion  term  in  the 
clarifier  model  equations  (Anderson  and  Edwards,  1981;  Lev  et  al.,  1986).  This  can  also 
give  a  realistic  solids  concentration  profile  that  can  be  used  as  the  basis  for  predicting 
sludge  blanket  height  (Hamilton  et  al.,  1992).  Recently,  the  Hamilton  et  al.  (1992) 
dispersion  model  was  modified  by  incorporating  dependence  of  the  dispersion  coefficient 
on  both  feed  velocity  and  local  solids  concentration  (see  Chapter  2). 

Work  reported  in  this  paper  improves  the  Hamilton  et  al.  (1992)  algorithm  for 
predicting  sludge  blanket  heights.  It  subsequently  shows  that  the  clarifier  model  can  be 
calibrated  without  measurements  of  the  suspended  solids  profile  along  the  depth  of  the 
clarifier.  The  calibration  uses  as  measure  of  fit  a  weighted  combination  of  the  error 
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squared  in  blanket  height,  the  error  squared  in  effluent  suspended  solids  concentration, 
and  the  error  squared  in  the  return  activated  sludge  (RAS)  concentration.  The  model  was 
calibrated  using  a  set  of  data  collected  at  the  Kanapaha  Water  Reclamation  Facility 
(KWRF)  in  Gainesville,  Florida,  then  validated  against  two  other  KWRF  data  sets 
collected  during  different  time  periods. 

Description  of  Clarifier  Model 
Dispersion  Coefficient  Expression 

As  is  the  case  with  Hamilton  et  al.  (1992),  the  model  incorporates  a  dispersion 
coefficient  D.  In  this  case,  however,  D  is  dependent  on  the  suspended  solids 
concentration.  The  dispersion  coefficient  expression  has  the  following  characteristics: 

a)  Decreasing  value  with  increasing  concentration,  at  high  solids  concentrations.  This 
could  perhaps  be  attributed  to  the  concomitant  increase  in  viscosity.  As  the  dispersion 
coefficient  should  always  be  positive,  an  exponential  decrease  towards  zero  is  reasonable. 

b)  A  constant  value  at  low  solids  concentrations. 

c)  A  smooth  transition  from  the  constant  dispersion  region  to  the  exponential  decay  region 
at  a  certain  concentration,  which  will  heretofore  be  referred  to  as  Ccrit. 

The  dependence  of  D  on  the  solids  concentration  C  is 

D(C)  =  Pmax[l+AC-Ccnt)]exp[-AC-Ccnt)]  forOC^ 

IA-  forC<Ccn,  U) 

Note  that  in  (1)  the  top  expression  yields  D(Ccrit)  =  D^  and  —  (Ccnt)  =  0,  i.e.,  we  have 
continuity  and  smoothness  at  C=Ccrit 


44 

In  the  model,  the  clarifier  is  divided  into  layers  numbered  from  top  to  bottom,  and 
the  dispersion  term  carries  material  from  a  layer  of  concentration  C;  to  a  layer  of 
concentration  Cm.  In  equation  (1)  the  concentration  C  that  affects  the  dispersion 
coefficient  is  taken  to  be  the  geometric  mean  of  Cj  and  Cm,  C,.i.,  .  It  should  be  remarked 

that  the  above  concentration-dependent  dispersion  coefficient  can  be  linked  to  the  Takacs 
et  al.  (1991)  model  as  shown  in  Chapter  2. 

The  dispersion  coefficient  was  also  correlated  to  velocity.  There  is  physical 
justification  for  this,  as  Taylor  (1953)  found  a  quadratic  dependence  of  the  dispersion 
coefficient  on  velocity.  It  is  shown  in  Chapter  2  that  fits  at  different  loadings  were 
markedly  improved  by  allowing  D™*  to  vary  as  a  function  of  the  feed  velocity  Vf.  The 
following  functional  dependence: 


f  Di+y(Vf-VM)2       ifVf>Vu 
D™  =  i  (2) 

I  Di  ifVf<Vf,i 


gave  good  results.  The  full  dispersion  coefficient  expression  is  given  by  equation  (1)  with 
Dmax  calculated  using  equation  (2). 

Settling  Velocity  Equation 

The  equation  proposed  by  Takacs  et  al.  (1991)  was  used  to  model  gravity  settling 
velocity  as  a  function  of  concentration.  This  equation  subtracts  an  exponential  term  from 
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the  commonly  used  Vesilind  (1968)  equation  resulting  in  decreasing  settling  velocity  with 
decreasing  concentration  in  the  low  concentration  region: 

Vs,,  =  min{V0  [expC-rXCj-Oun))  -  expC-bpCCi-Cs.))] ,  V™  }  (3) 

In  this  equation,  VS)i  is  the  gravity  settling  velocity  from  layer  i,  C;  is  the  layer  suspended 
solids  concentration,  Cm,,  is  the  nonsettleable  suspended  solids  concentration,  V0  and  b  are 
the  Vesilind  (1968)  settling  parameters,  bp  is  a  settling  parameter  characteristic  of  low 
suspended  solids  concentrations,  and  Vmax  is  the  highest  settling  velocity  achieved  by 
sludge  floes. 

Model  of  KWRF  Clarifier 

A  schematic  of  the  KWRF  test  clarifier  is  presented  in  Figure  3-1 .  The  model 
accounts  for  the  presence  of  an  annular  baffle  (shroud)  in  the  upper  cylindrical  section  of 
the  clarifier  (from  0  to  2.44  m  below  the  water  surface).  The  cross-sectional  area 
available  for  overflow  (Ac)  in  this  region  is  72%  of  the  cross-section  area  (A)  in  the  lower 
cylindrical  section  (between  the  shroud  and  the  cone)  of  the  clarifier. 

The  bottom  of  the  clarifier  consists  of  a  conical  section  from  which  RAS  is  removed 
by  a  hydraulic  suction  system.  Because  the  pipe  intakes  are  spaced  at  different  depths 
along  the  sloping  bottom,  the  concentrations  of  sludge  withdrawn  at  the  different  depths 
will  be  substantially  different.  Therefore,  sludge  removal  should  not  be  modeled  as  if  it 
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Fig.  3-1.  Schematic  diagram  of  clarifier  at  KWRF. 
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were  withdrawn  from  a  single  layer.  As  Figure  3-1  shows,  the  conical  section  is  divided 
into  p  layers  and  the  cylindrical  section  into  n  layers.  In  the  model,  sludge  is  removed 
from  each  layer  of  the  conical  section  at  a  rate  equal  to  the  change  in  cross-sectional  area 
from  the  top  to  the  bottom  of  the  layer  multiplied  by  the  underflow  velocity: 

Qr,i  =  (Ai-Ai+1)Ub  (4) 

Here  Qr>j  is  the  rate  of  flow  removed  from  layer  i,  A;  is  the  area  at  the  top  of  a  layer  in  the 
conical  section,  and  Ai+i  is  the  area  at  the  bottom  of  the  layer.  The  underflow  velocity  Ub 
in  the  lower  sections  of  the  clarifier  is  the  ratio  of  the  sum  of  the  RAS  and  waste  activated 
sludge  flows,  Qu,  to  the  area  of  the  lower  cylindrical  section,  A.  Note  that 

n+p 

Qu  =  Z  Q     and  that  allocation  of  the  withdrawal  flows  according  to  eq.  (4)  maintains  a 

i=n+l 

constant  underflow  velocity  throughout  the  conical  section.  The  RAS  concentration  will 
be  the  flow-weighted  average  of  the  concentrations  of  layers  in  the  conical  section 

Let  the  thickness  of  the  layers  be  Sz  and  let  the  layers  be  numbered  from  top  to 
bottom.  Mass  balances  on  suspended  solids  give  the  equations  below 

For  the  top  layer  (i  =  1): 

5z  dCi/dt  =  Uq  C2  -  UqCi  -  VoC,  +  D,,2  (C2  -  C,)/  5z  (5) 

For  the  \th  layer  in  the  region  above  the  feed  layer: 

5z  dCi/dt  =  Uq  Chi  -  UqCi  +  Vs,i.,Ci.i  -  Vs,iC,  -  Dmj  (Ci-Ci.i)/5z  +  D;,i+1  (Ci+1-Ci)/5z     (6) 
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For  the  layer  receiving  clarifier  feed: 

6z  dCi/dt  -  (Uq+Ub)  Cf  -  Uq  C,  -  UbCi  +  V^C,.,  -  V^Q  -  Di.u  (Ci-Cj.i)/8z 

+  Di,i+1  (Ci+r  Ci)/6z  (7) 

For  the  \th  layer  in  the  region  below  the  feed  layer: 

8z  dCi/dt  -  Ub  Ci-i  -  Ub  Q  +  Vs.-iCi.,  -  V^Q  -  Df.u  (CrQ.O/Sz  +  Diii+i  (Ci+1-Ci)/8z   (8) 

And  finally  for  the  bottom  layer: 

5z  dCi/dt  =  Ub  Ci.i  -  UbCi  +  Vs,,.iCi.i  -  Dm,;  (Q-Ci.,)/5z  (9) 

Here  Uq  is  the  overflow  velocity  (=  Qe/Ac  in  the  upper  cylindrical  region  or  QJA  in  the 
lower  cylindrical  section  with  Qe  being  the  effluent  flowrate),  Ub  the  underflow  velocity 

(-  Qu/Ac  in  the  upper  cylindrical  section  or  QJA  below  the  shroud),  and  Dj_i,i  =  D(Ci.ij). 
To  calculate  the  location  of  the  feed  layer,  an  initial  position  for  the  feed  layer  was 
assumed  and  the  concentration  profile  calculated.  The  location  of  the  feed  layer  was  then 
chosen  as  the  layer  above  the  first  layer  having  a  concentration  greater  than  the  feed 
concentration.  This  is  consistent  with  observations  of  density  current  flows  in  prototype 
scale  clarifiers  (Andersen,  1945).  Subsequently,  the  concentration  profile  was 
recomputed,  and  the  feed  layer  reassigned.  The  latter  steps  were  repeated  until 
convergence  was  achieved.  In  rare  instances  when  the  above  procedure  did  not  converge 
to  a  single  layer  and  instead  began  to  oscillate  between  two  adjacent  layers,  the  higher  of 
the  two  layers  was  chosen  as  the  feed  layer. 
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Sludge  Blanket  Algorithm 

The  sludge  blanket  algorithm  was  modified  from  Hamilton  et  al.  (1992),  who 
calculated  the  blanket  height  as  the  height  corresponding  to  the  maximum  rate  of  change 
in  the  slope  Pi  of  the  solids  concentration  versus  depth  profile.  The  modified  algorithm 
uses  in  place  of  Pj  a  relative  concentration  slope: 

(C,  +  C,_,)/2 
which  is  the  concentration  slope  divided  by  the  average  concentration  between  the 
adjacent  layers.  This  gives  higher  predicted  blanket  heights  than  those  calculated  with  the 
Hamilton  et  al.  (1992)  algorithm. 

The  height  hi  (measured  from  the  bottom  of  the  conical  section)  of  the  interface 
between  model  layers  having  the  maximum  relative  concentration  slope  R;  is  located,  then 
a  quadratic  interpolating  polynomial  is  used  to  find  a  smooth  curve  passing  through  this 
point  and  two  adjacent  points  (R.i,  hj.i;Ri+i,  hi+i).  When  the  maximum  Rj  is  at  the 
interface  between  the  bottom  layer  and  the  layer  above  it,  R+i  (at  h;+1  =  0)  cannot  be 
calculated  by  eq.  (10).    It  is  then  assigned  a  value  of  zero  to  ensure  that  a  positive  blanket 
height  is  obtained.  The  sludge  blanket  height  (SBH)  is  calculated  as  the  location  of  the 
maximum  of  the  interpolated  polynomial.  This  results  in 

SDII_~(m2'R-1+  m22Rl  +  m"R»'>  (ii) 

2(m11Rll+m12R1+ml3Rl+1) 

where  m,j  is  the  ij*  element  of  the  matrix 
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hM    Ki    ] 


-1 


M  =      h?       h,      1  (12) 


h?+i     h,+i     1 


Materials  and  Methods 
Clarifier  Loading  Tests 

Clarifier  loading  tests  were  conducted  at  the  38,000  m3/d  (10  Mgal/d)  Kanapaha 
Water  Reclamation  Facility  (KWRF)  using  one  of  the  plant's  four  secondary  clarifiers. 
The  clarifier  was  28.96  m  (95  ft)  in  diameter  with  a  3.66  m  (12  ft)  sidewall  depth  and  a  4.6 
m  (15  ft)  depth  at  the  center. 

During  a  loading  test,  influent  and  RAS  flow  rates  were  controlled  at  selected  values. 
Sludge  blanket  height  was  measured  every  1 5  min  using  a  5  cm  ID  transparent  plastic  tube 
(Sludge  Judge,  NASCO  Inc.,  Ft.  Atkinson,  Wisconsin).  Data  collection  was  continued  for 
two  hours  after  reaching  a  steady  blanket  level.  Tests  where  a  continuously  rising  blanket 
was  observed  were  ended  when  the  sludge  blanket  approached  the  effluent  weirs.  Influent 
mixed  liquor  was  sampled  hourly  throughout  each  test,  whereas  clarifier  effluent,  RAS, 
and  waste  activated  sludge  (WAS)  were  sampled  hourly  once  a  steady  blanket  level  was 
reached  or  when  the  sludge  blanket  approached  the  effluent  weirs.  Samples  were  stored 
on  ice  until  analysis  for  suspended  solids  later  in  the  day. 

To  determine  concentration  profiles,  samples  were  collected  at  0.61  m  intervals 
through  the  clarifier  depth.  The  sampling  location  was  halfway  between  the  shroud  and 
the  peripheral  wall  of  the  clarifier.  An  additional  location  was  inside  the  shroud. 
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Batch  Settling  Test  Procedure 

Batch  sludge  settling  tests  were  conducted  in  a  water-bath  enclosed  six-column 
settling  apparatus  (Wahlberg  and  Keinath,  1986).  RAS,  mixed  liquor,  and  clarifier  effluent 
were  mixed  in  selected  ratios  to  obtain  suspended  solids  concentrations  ranging  from  2000 
to  14000  mg/L.  Columns  were  mixed  for  5  min  prior  to  each  test  using  compressed  air. 
Interface  height  in  each  column  was  recorded  every  two  minutes  until  the  compression 
phase  was  reached. 

Settling  velocity  (Vs)  at  each  initial  suspended  solids  concentration  (C)  was  found  by 
least  squares  linear  regression  on  the  initial  linear  portion  of  the  interface  height  versus 
time  curve.  Settling  parameters  V0  and  b  in  the  expression  (Vs  =  V0  e*bc  )  (Vesilind,  1968) 
were  found  by  least  squares  linear  regression  on  the  logarithms  of  the  settling  velocities 
and  the  corresponding  sludge  concentrations  (Daigger,  1995).  Four  settling  tests  were 
carried  out  during  experimental  period  A,  six  settling  tests  during  period  B,  and  seven 
during  period  C.  (A  settling  test  refers  to  the  six  settling  trials  carried  out  simultaneously 
in  the  multi-column  apparatus.)  Settling  test  data  were  screened  by  comparing  the  V0  and 
b  from  each  individual  test  to  the  preliminary  mean  V0  and  b  for  all  tests  in  the  respective 
test  period.  Tests  with  Vo  or  b  more  than  2.0  standard  deviations  away  from  the 
preliminary  mean  parameter  values  from  that  period  were  rejected.  Accepted  data  from 
the  respective  test  periods  (only  one  test  was  rejected)  were  pooled  and  parameters  of  the 
Vesilind  expression  were  found  by  linear  regression.  Parameters  were  V0  =  149.9  m/d  and 
b  =  3.921  x  10"4  L/mg  for  period  A  V0  =  152.5  m/d  and  b  =  3.213  x  10"4  L/mg  for  period 
B,  and  V0  =  182.9  m/d  and  b  =  3.055  x  10"4  L/mg  for    period  C. 
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Additional  Measurements 

The  mean  of  grab  samples  from  the  clarifier  effluent  (collected  for  use  in  the  batch 
settling  tests)  was  used  as  the  nonsettleable  suspended  solids  concentrations,  Cmin,  for 
simulations.  These  values  were  3.9  mg/L  for  period  A  and  2.4  mg/L  for  both  periods  B 
and  C.  Samples  for  total  suspended  solids  analysis  were  filtered  through  glass  fibre  filters 
having  an  average  pore  size  of  1.2  pm  (Whatman  GF/C).  Filtered  residues  were  dried  to 
constant  weight. 

Parameter  Estimation 

Model  parameters  (Vmax,  bp,  Dh  y,  Vf,i)  were  estimated  using  the  Levenberg- 
Marquardt  algorithm  (Marquardt,  1963;  Press  et  al.,  1989)  with  scaling  according  to 
Cuthbert  (1987).  The  objective  function  weighted  three  components:  the  sum  of  the 
squares  of  the  relative  errors  in  sludge  blanket  height,  the  sum  of  the  squares  of  the 
relative  errors  in  RAS  concentration,  and  the  sum  of  the  squares  of  the  relative  errors  in 
effluent  concentration.  The  three  components  were  weighted  80%:  10%:  10%, 
respectively,  since  the  blanket  height  provides  the  greatest  amount  of  information  about 
the  shape  of  the  concentration  profile. 
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Results  and  Discussion 
Clarifier  Loading  Tests 

A  total  of  43  clarifier  loading  tests  were  carried  out  during  three  experimental 
periods  (A,  B,  C)  over  a  time  span  of  6  months.  Out  of  13  tests  in  period  A,  nine  reached 
steady-state  conditions  as  judged  from  steady  blanket  levels.  Eleven  of  the  15  tests  in 
period  B  and  10  of  the  15  tests  in  period  C  also  reached  steady-state  conditions.  Mass 
balance  closure  errors  under  steady-state  conditions  ranged  from  0.9%  to  18.6%,  with  a 
median  error  of  9.0%.  Overall  mean  influent  suspended  solids  concentrations  ranged 
between  3200  and  4170  mg/L,  whereas  effluent  suspended  solids  concentrations  under 
steady-state  conditions  ranged  from  1.6  mg/L  to  8.7  mg/L.  Effluent  suspended  solids 
were  not  correlated  with  total  solids  loading  (P<0.05).  Steady-state  RAS  concentrations 
ranged  between  6330  and  14600  mg/L. 

Model  Calibration 

The  clarifier  model  was  calibrated  using  steady-state  data  (sludge  blanket  height, 
effluent  suspended  solids  and  RAS  concentrations)  from  test  period  C.  The  estimated 
model  parameters  were  Vmax  =  172.0  m/d,  bp  =  2.7587  x  10"2  L/mg,  Di  =  4.835  m2/d, 
y  =  2.500  x  10"2  d,  and  Vf,i  =  32.88  m/d.  Model  profiles  exhibited  a  zone  of  rapidly 
changing  suspended  solids  concentration  between  regions  of  slowly  changing 
concentration  near  the  top  and  bottom  of  the  clarifier,  and  were  generally  consistent  with 
measured  data  (Fig.  3-2).  As  Table  3-1  shows,  the  model  tended  to  overestimate  RAS 
concentration  slightly  (median  error  =  4.7%).  The  error  in  RAS  concentration  predictions 
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is  essentially  set  by  the  mass  balance  closures  of  the  respective  experiments.  The  model 
tended  to  underestimate  effluent  suspended  solids  concentration  by  as  much  as  3.3  mg/L, 
and  the  correlation  of  the  model  estimates  to  the  measurements  was  poor.  This  is  not 
surprising  since,  in  the  tests,  effluent  suspended  solids  were  not  correlated  with  clarifier 
loading.  (This  holds  since  experiments  were  never  carried  out  to  loss  of  blanket.) 

Sludge  blanket  heights  are  represented  in  Figure  3-2  by  vertical  lines,  light  solid  for 
model-calculated  heights  and  heavy  solid  for  measured  heights.  Model  sludge  blanket 
heights  fall  in  the  region  of  rapidly  changing  concentration,  as  indicated  by  the  model 
profile.  Model  blanket  heights  differed  from  measurements  by  no  more  than  0.21  m,  and 
were  within  0.08  m  of  measurements  in  most  tests.  A  linear  regression  of  model- 
calculated  sludge  blanket  heights  to  the  measurements  gave  a  slope  of  1.04  (Fig.  3-3). 

Model  Validation 

The  model  as  calibrated  using  test  period  C  data,  along  with  V0,  b,  and  Cmin  measured 
in  each  period,  was  run  to  predict  sludge  blanket  heights  of  tests  in  periods  A  and  B.  The 
relationship  between  predicted  and  measured  heights,  for  the  tests  in  which  a  steady 
blanket  was  observed,  is  shown  in  Figure  4.  In  period  A,  for  all  nine  tests  the  model 
blanket  heights  within  0.43  m  (Table  3-2)  with  median  absolute  error  of  0. 16  m.  For  ten 
of  the  eleven  tests  of  period  B,  the  model  predicted  blanket  height  within  0.36  m  of  the 
measurement.  The  one  exception  was  the  test  with  the  lowest  underflow  velocity  (case  B- 
7),  in  which  the  model  overpredicted  the  height  by  2. 1  m.  The  median  absolute 
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Fig.  3-3.  Comparison  of  predicted  blanket  heights  to  measured  values  from  test  period  C. 
(Model  was  calibrated  using  blanket  height  data  in  addition  to  RAS  and  effluent 
suspended  solids  concentrations  from  test  period  C.) 
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Fig.  3-4.  Comparison  of  predicted  blanket  heights  to  measured  values  from  test  periods 
A  and  B.  (Model  parameters  from  calibration  to  test  period  C  data  were  used  in  the 
simulations.) 
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Table  3-2.  Comparison  of  predicted  effluent  concentrations,  RAS  concentrations, 
and  sludge  blanket  heights  to  measured  values  for  test  periods  A  and  B 


Measurements 

Predictions  * 

effluent 

RAS 

blanket 

effluent 

RAS 

blanket 

Case 

cone. 

cone. 

height 

cone. 

cone. 

height 

(mg/L) 

(mg/L) 

(m) 

(mg/L) 

(mg/L) 

(m) 

A-l 

7.6 

11358 

1.86 

4.5 

11589 

1.95 

A-2 

7.2 

6691 

1.35 

4.0 

7347 

1.19 

A-3 

5.4 

8313 

1.41 

4.3 

8872 

1.39 

A-4 

6.9 

8313 

1.97 

4.4 

9138 

1.54 

A-5 

7.2 

7152 

1.56 

4.3 

7859 

1.29 

A-6 

8.7 

6332 

1.59 

4.3 

7404 

1.29 

A-7 

6.9 

9790 

2.29 

5.1 

10219 

2.18 

A-8 

5.9 

9075 

2.01 

4.9 

9692 

1.86 

A-9 

3.0 

8806 

2.16 

4.9 

9525 

1.79 

B-l 

3.6 

8383 

1.11 

3.2 

7830 

1.15 

B-2 

3.9 

9155 

1.52 

3.8 

10376 

1.49 

B-3 

3.6 

9679 

1.98 

4.8 

11186 

2.00 

B-4 

3.2 

9681 

2.36 

4.8 

11165 

2.00 

B-5 

1.9 

10170 

2.23 

4.8 

10975 

1.99 

B-6 

2.7 

9271 

1.21 

3.2 

10152 

1.25 

B-7 

3.6 

12788 

2.22 

40.8 

15340 

4.34 

B-8 

1.9 

7774 

2.08 

5.4 

9216 

2.02 

B-9 

1.6 

7324 

1.89 

5.0 

7721 

1.77 

B-10 

2.5 

9040 

2.23 

5.4 

9053 

2.01 

B-ll 

4.9 

7470 

1.20 

3.2 

8588 

1.19 

*Model  parameters,  with  the  exception  of  V0,  b,  and  Cmi„,  were  those  found  from 
test  period  C  fit 
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error  of  the  eleven  tests  was  0.06  m.  The  linear  regression  slope  for  all  the  predictions  in 
periods  A  and  B  was  1.11  (Fig.  3-4). 

Table  3-2  also  shows  the  predictions  of  RAS  and  effluent  suspended  solids 
concentrations  for  the  tests  in  periods  A  and  B.  The  model  tended  to  over  predict  RAS 
concentrations  in  a  manner  consistent  with  the  mass  balance  closures  of  the  respective 
experiments.  The  range  of  predicted  effluent  suspended  solids  concentrations  for  periods 
A  and  B,  excluding  case  B-7,  was  3.2  -  5.4  mg/L  with  a  median  (including  B-7)  of  4.65. 
In  comparison,  the  range  of  the  measured  suspended  solids  concentrations  was  1.6-8.7 
mg/L  with  median  3.75  mg/L.  Other  than  having  values  of  the  correct  order  of 
magnitude,  the  model  predictions  do  not  correlate  well  with  the  measured  values. 

Tables  3-3a,  3-3b,  and  3-3c  give  the  results  of  the  clarifier  tests  and  classify  each  test 
according  to  whether  the  applied  loading  was  acceptable  or  not  under  the  test  operating 
conditions.  Since  field  tests  could  not  be  carried  out  to  the  point  of  actual  failure,  the  end 
result  was  judged  on  the  basis  of  sludge  blanket  behavior.  Experimentally  observed  sludge 
blankets  that  were  rising  at  a  rapid  rate  (>  0. 1  m/h)  near  test  termination  were  considered 
to  be  the  result  of  clarifier  overloading,  whereas  steady  or  falling  blanket  levels  near  test 
termination  were  considered  to  result  from  acceptable  loading  levels.  Tests  ending  in 
slowly  rising  blankets  were  considered  inconclusive. 

Also  reported  in  Table  3-3  are  the  results  of  simulations  with  the  present  model 
under  the  loading  and  operating  conditions  of  each  test  run.  V0,  b,  and  Cmin  were  the 
measured  values  for  each  of  the  three  test  periods,  whereas  all  other  model  parameters 
were  those  found  in  test  period  C.  Model  results  were  classified  according  to  the 
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predicted  effluent  suspended  solids.  Effluent  concentrations  of  less  than  20  mg/L  were 
considered  to  represent  acceptable  loading  levels.  "Overloaded"  cases  were  characterized 
by  effluent  concentrations  in  excess  of  70  mg/L.  Predictions  of  the  Takacs  et  al.  (1991) 
model,  which  was  considered  to  be  the  best  available  by  Grijspeedt  et  al.  (1995),  as  well  as 
limiting  solids  flux  theory  (Coe  and  Clevenger  1916;  Yoshioka  et  al.  1957)  are  given  in  the 
table  for  comparison.  The  sludge  blanket  algorithm  was  incorporated  in  a  program 
implementing  the  Takacs  et  al.  (1991)  model,  and  the  model  was  calibrated  to  KWRF  test 
period  C  data  using  procedures  identical  to  those  described  above.  Simulations  of  the 
other  test  periods  were  performed  with  the  test  period  C  fitted  parameters  and  V0,  b,  and 
Cmin  for  the  given  period.  The  total  limiting  solids  flux  was  calculated  based  on  the 
measured  V0  and  b  values  for  the  appropriate  test  periods. 

Out  of  40  conclusive  loading  tests,  the  present  model  correctly  predicted  the 
outcome  of  37.  Two  experimental  successes  were  predicted  as  failures  (cases  A-l  1  and 
B-7),  whereas  one  experimental  failure  was  predicted  as  a  success  (case  C-15),  but  with  a 
relatively  high  blanket  height  and  effluent  solids  concentration.  In  contrast,  the  Takacs  et 
al.  (1991)  model  predicted  success  for  all  cases  (which  is  also  consistent  with  limiting 
solids  flux  theory),  and  thus  incorrectly  predicted  the  eight  experimental  failures. 

Conclusions 
The  model  utilized  in  the  present  work  introduces  an  algorithm  for  locating  the  top 
of  the  sludge  blanket  based  on  the  point  of  greatest  relative  concentration  slope.  This 
algorithm  is  computationally  efficient  and  reliable  in  matching  experimentally  measured 
blanket  heights.  Incorporation  of  the  blanket  algorithm  in  the  model  enables  calibration 


65 

using  measured  blanket  heights  instead  of  concentration  profiles  throughout  the  clarifier 
depth.  The  model  validation  carried  out  in  this  research  was  based  on  extensive  full  scale 
plant  data  sets.  Based  on  these  data,  the  model  developed  in  this  work  appears  to  be  more 
reliable  than  limiting  solids  flux  theory  or  the  model  of  Takacs  et  al.  (1991)  in  predicting 
clarifier  failure  due  to  solids  overloading. 


CHAPTER  4 
CONCLUSIONS 


As  a  result  of  the  modeling  and  experimental  work  carried  out  in  this  study,  the 
following  conclusions  can  be  drawn: 

•  Clarifier  models  incorporating  a  constraint  on  gravity  flux  have  been  shown  to  provide 
good  fits  to  experimental  concentration  profiles,  but  the  flux  constraint  is  dependent 
upon  the  level  of  model  discretization.  The  flux  constraint  therefore  effectively 
disappears  as  the  level  of  discretization  is  increased. 

•  The  gravity  flux  constraint  can  be  recast  as  a  concentration-dependent  dispersion  term 
which  improves  the  model's  ability  to  fit  experimental  data  as  the  level  of  model 
discretization  is  increased. 

•  The  ability  of  the  model  to  fit  concentration  profiles  collected  in  steady-state  full-scale 
clarifier  loading  tests  over  a  range  of  solids  loading  was  substantially  improved  by 
inclusion  of  dependence  on  influent  velocity  in  the  dispersion  function. 

•  Incorporation  of  the  algorithm  for  the  determination  of  sludge  blanket  height  in  the 
model  enables  calibration  using  measured  blanket  heights  instead  of  concentration 
profiles. 

•  Model  validation  carried  out  in  this  research  was  based  on  extensive  full-scale  plant 
data  sets.  The  model  was  calibrated  using  data  from  nine  steady-state  clarifier  loading 
tests  from  one  experimental  period,  and  was  validated  against  data  from  twenty 
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steady-state  clarifier  loading  tests  from  two  different  experimental  periods.  The  model 
was  further  tested  in  simulations  of  fourteen  additional  loading  tests  from  all  three 
experimental  periods. 

•  Model  predictions  of  sludge  blanket  heights  and  underflow  suspended  solids 
concentrations  were  generally  good.  Effluent  suspended  solids  concentration 
predictions  were  of  the  correct  order  of  magnitude  but  did  not  correlate  with  measured 
concentrations. 

•  The  model  developed  as  a  part  of  this  research  outperformed  both  total  limiting  solids 
flux  theory  and  the  gravity-flux-constraining  model  in  the  prediction  of  clarifier  failure 
due  to  solids  overloading. 
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